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Effects of lattice distortion on the magnetic ground state of YTiOs and LaiOs are investigated on 
the basis accurate tight-binding parametrization of the tig electronic structure extracted from the 
local-density approximation. The complexity of these compounds is related with the fact that the 
£20-level splitting, caused by lattice distortions, is comparable with the energies of superexchange 
and spin-orbit interactions. Therefore, all these interactions are equally important and should be 
treated on an equal footing. The Hartree-Fock approximation fails to provide a coherent description 
simultaneously for YTiOs and LaTiOs, and it is essential to go beyond. 
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Among the large variety of transition-metal perovskite 
oxides, YTi0 3 (YTO) and LaTi0 3 (LTO) have received 
a particular attention. Both are regarded as prototypical 
examples of Mott-Hubbard insulators. It appears, how- 
ever, that these, formally isoelectronic compounds (hav- 
ing one 3d electron in the triply-degenerate t2 g shell), 
exhibit very different magnetic properties: YTO is a fer- 
romagnet, whereas LTO is a three-dimensional (G-type) 
antiferromagnet. Another puzzling feature is the nearly 
isotropic magnon spectrum, observed both in YTO and 
LTO despite a noticeable orthorhombic distortions^ 

Owing to the fractional population of the t^ g mani- 
fold, the orbital degrees of freedom are expected to play 
a very important role and affect the magnetic properties. 
However, the theories proposed in this context crucially 
depend on several factors, and there are two points of 
view which are currently discussed in the literature, (i) 
The first one is based on the generalization of the su- 
perexchange (SE) theory of spin and orbital interactions 
between degenerate ti g levels. It starts with the spin- 
orbital SE model by Kugel and Khomskii (K&K), 3 and 
exploits the idea of orbital fluctuations, which are inher- 
ent to this model. 4 (ii) The spin and orbital structure 
is fully determined by lattice distortions, which lift the 
orbital degeneracy The role of (relativistic) spin- 
orbit (SO) interaction has been also emphasized. 10 

Therefore, there are two important questions, which 
can be clarified on the basis of electronic structure calcu- 
lations, (i) What is the effect of the lattice distortion on 
the electronic structure of YTO and LTO? Particularly, 
how do the t<i g levels split by this distortion? (ii) What 
is the hierarchy between the ^-level splitting, the SE 
interaction energy, and the SO coupling? 

The SE interaction in the bond i-j is basically the ki- 
netic energy gain, which is acquired by the t^g electron 
occupying the atomic orbital |i) at the site i in the pro- 
cess of virtual hoppings into the subspace of unoccupied 
orbitals Pj at the (neighboring) site j, and vice versa>Ai 



(ifaPjtjili) + (i<-j) 
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where a and a' are the spin states associated with the 



sites i and j, respectively, and the transfer interactions 
tij are allowed only between orbitals with the same spin. 
For the nearest-neighbor interactions in the perovskite 
lattice, it is sufficient to consider two collinear configu- 
rations, gg'—\\ and ||, and select the ones which min- 
imize the total energy gain £t=^ Xaj £ ij a • ^ n the case 
of the antiferromagnetic (AFM) alignment, crcr'=]l, all 
orbitals with the spin j at the site j are located in the 
unoccupied part of the spectrum and available for the 
hoppings. Therefore, ij=l and a^ = (i|tytji|i) + (i<-»j). 
In the ferromagnetic (FM or F) case, <j<j'=^, the oc- 
cupied orbital |j) should be excluded from the subspace 
Pj. This yields Pj=l— |j)(j| and =a]^ — Aay, where 

Aaij=2 |(i|tij|j)| 2 . AE aa ' is the on-site Coulomb inter- 
action between two 3d electrons, which also depends on 
the spin state: AE n =U while AP TT =/7-J, where U 
is the Coulomb repulsion and J is the intra-atomic ex- 
change coupling.— Because of this J, the "orthogonal" 
orbitals, which do not interact via the kinetic energy 
term, (i|tij|j)=0, tend to stabilize the FM structure. In 
the opposite limit Aa^a^ , the FM alignment does not 
lead to any energy gain, and the coupling will be AFM. 

The alternation of occupied orbitals at different atomic 
sites (the orbital ordering - 00) should be found varia- 
tionally and minimize St- This is the basic idea of the 
K&K theory. 3 The orbital interactions have the same ori- 
gin as the spin SE. Therefore, the energy gain associated 
with the 00 is of the order of Et~1/U, and there is 
a strong interplay between spin and orbital degrees of 
freedom. In the degenerate case, one can always find 
some orthogonal configuration of the occupied orbitals, 
which in the single-determinant Hartree-Fock (HF) ap- 
proach corresponds to the FM ground state (GS). How- 
ever, the HF solutions remain degenerate with respect to 
some number of orbital configurations. This degeneracy 
leaves a room for orbital fluctuations, which may alter 
the HF conclusion about the form of the magnetic GS. 4 

An alternative mechanism of the 00 is the lattice dis- 
tortions, which lifts the orbital degeneracy and acts as 
an external field constraining the form of occupied or- 
bitals in Eq. JQ. Since the orbital degeneracy is lifted, 
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the HF approach may be justified. 6,8 This mechanism is 
proportional to the electron-phonon coupling, and will 
dominate over the K&K SE mechanism in the large- U 
limit. Then, the 00 does not depend on the magnetic 
state and the mapping onto the Heisenberg model yields 
the following expression for J^=^sj^— sj^): 



all J/U - Aau/ ajj- 
2 U — J 



Jy = 



U 



(2) 



which can be both FM and AFM, depending on the ratio 
of J/U and Aay/oJp 

Let us consider the second scenario and assume that 
all relevant interactions can be described in the ba- 
sis of some local t2 g orbitals \X), \Y), and asso- 
ciated with the Ti sites. Then, the occupied orbital 
at the site 1 (see Fig. [Q can be searched in the form 
|l)=sin#cos(/>|X)-hsin#sin0|y)+cos#|Z), and the ones 
at the sites 2 and 3 are automatically generated from 
|1) using the symmetry operations of the group (the 
180° rotations around the orthorhombic a and c axes, 
respectively). In principle, 6 and (j) are uniquely deter- 
mined by the lattice distortion. 

However, it is sometimes tempted to approach the 
problem from the opposite side, 7 and find 9 and (j) 
from the condition Ji2=Ji3, suggested by recent neu- 
tron scattering studies m£ In order to illustrate this idea, 
let us consider a simplified model and choose \Y), 
and \Z) as \yz), \zx), and \xy), respectively, in the 
cubic coordinate frame shown in Fig. ^ The trans- 
fer interactions are parameterized according to Slater 
and Koster (S&K): i.e., the only nonvanishing ma- 
trix elements along z are tf 3 x =t\^ '=t, etc. Then, 
it is easy to verify that the condition Ji2=Ji3 leads 
to the following OO: 13 \l)=\2)=-±(\xy)+\yz)+\zx)) : 

\3) = \4) = -j=(\xy) — \yz) — \zx)), which does not depend on 
J/U. This is precisely the 00 proposed in Ref. 0. It 
is compatible with the orthorhombic D^® symmetry, and 
corresponds to some local trigonal distortion, caused by 
either oxygen or La displacements. 8 

This result, however, prompts several new questions, 
(i) The magnetic coupling is expected to be AFM for all 
reasonable values of J/U. Therefore, this would explain 
the experimental situation in LTO, but not in YTO. (ii) 
It is not clear whether this 00 is compatible with the ac- 
tual experimental distortion observed in LTO. Note that 
in the D^® group, only inversion centers coincide with 
the Ti sites. Therefore, the local ^-level splitting is con- 
trolled by 5 independent parameters, which may include 
both trigonal and Jahn- Teller modes. All distortions are 
formally equivalent, at least from the viewpoint of 
symmetry, and a priori there is no reason why the par- 
ticular trigonal mode should dominate. In addition to 
the ^2^-level splitting, the crystal distortion may also af- 
fect the transfer interactions through the buckling of the 
Ti-O-Ti bonds. 14 (hi) What are the roles of the K&K 
mechanism and the SO interaction? Are they totally 
quenched by the lattice distortion, as it was suggested 





FIG. 1: t2^-electron densities obtained in Hartree-Fock cal- 
culations after including the spin-orbit interaction, x, y, and 
z are the cubic axes, a, b, and c are the orthorhombic axes. 



in Refs. The situation should be carefully checked, 
and it is important to turn to first-principles calculations, 
which automatically include all these ingredients. 

We use the linear-muffln-tin-orbital (LMTO) 
method^ and employ the tight-binding (TB) 
parametrization of the t^g bands, obtained in the 
local-density approximation (LDA) for the experimental 
crystal structures^ The latter step is achieved through 
the downfolding procedure. A similar analysis has been 
undertaken in Ref. (i) Each LMTO eigenvector is 
divided in two parts: which is expanded over the 
local t2 g orbitals \X), \Y), and \Z) at each Ti site, 
and |r), which is expanded over the rest of the basis 
functions. The corresponding secular equation, which 
holds for the LMTO Hamiltonian 7Y, is given by 



(Htt-E)\t) 

Hrt\t) + (H 



H tr \r) 
-E)\r) 



0, 
0. 
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(ii) By eliminating \r) from Eq. (0J one ob- 
tains an effective independent Hamiltonian: 
nf(E)=n tt -n tr (n rr -E)- 1 n rt , where \t) obeys the 
condition (t\S\t)=l and S(E)=l+H tr (H rr -E)- 2 H rt . 
(hi) The TB parameters t=||ty|| are obtained after the 
orthonormalization of the vectors \t)^\t)=S 1 ^ 2 \t): 



t\E) = S- 1/2 (E)H e *(E)S- i/2 (E), 



eff / 



-1/2/ 



(5) 



Finally, E is fixed to the center of the t<i g band. 

The choice of the local t<i g orbitals is somewhat am- 
biguous. In our case we first calculated the site-diagonal 
elements of the density matrix in the basis of all Ti(3d) 
orbitals and taking into account the contributions of only 
the t2 g bands shown in Fig. [3 This yields the 5x5 matri- 
ces at each Ti site. Then, we assign three most populated 
orbitals obtained after the diagonalization of these ma- 
trices to \X), \Y), and \Z). 

The mapping onto the TB model is nearly perfect and 
well reproduces the behavior of LMTO bands (Fig. [5J). 
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FIG. 2: Energy bands obtained in LMTO calculations and 
after tight-binding (TB) parametrization. 



Then, the site-diagonal elements of t\j describe the 
crystal-field (CF) splitting caused by lattice distortions, 
and the off-diagonal elements have a meaning of trans- 
fer interactions. Thus, we are ready to calculate the SE 
interactions in the strong coupling limit, assuming that 
the form of occupied orbitals is solely determined by 
and using these orbitals in subsequent calculations of aj^ 
and Aaij. The results are summarized in Table Q] 

(i) The CF splitting is larger in YTO, mainly due to 
the Jahn- Teller distortion, 9 which is reflected in the up- 
ward shift of one of the t2 g levels™ The CF splitting in 
LTO is not particularly strong (in fact it is considerably 

weaker than the model estimates presented in Refs. Isfci)- 

1 1 

The inter-atomic interactions aJj are larger in the less 
distorted LTO, that well correlates with the larger t^g 
bandwidth (Pig. EJ). (ii) Both compounds exhibit certain 
tendency to A-type antiferromagnetism, which is espe- 
cially strong in YTO: since Aai2~0 and Aa^^a^, the 
bonds 1-2 and 1-3 are expected to be FM and AFM, 
respectively, for all physical values of J/U. Therefore, 
the crystal distortion alone cannot explain the FM GS 
of YTO. 18 The situation is somewhat milder in LTO 
where the experimental G-type AFM ordering can be 
stabilized for J ' jJJ "<Aai2/ 'a|2~0.37 '. 12 However, even in 
this case the interatomic magnetic interactions are ex- 
pected to be anisotropic, (iii) Realistic estimates for the 
on-site Coulomb interaction U in the t^g band typically 
vary from 3.2 eV, suggested by constraint-LDA calcula- 
tions and taking into account the empirical screening by 
the e g electrons^ to 4.4 eV suggested by photoemission 
studies™*^ The intra-atomic exchange coupling can be 
estimated as J-0.9 eV. 1Q i 19 i 2Q Therefore, £t can be as 
large as 10-40 meV per one Ti site. This value can be 
used as a rough estimate for the 00 stabilization energy 
caused by SE interactions, which is comparable with the 
CF splitting. Therefore, the K&K mechanism remains 



TABLE I: The crystal-field (CF) splitting of the t 2g states 
(in meV) and parameters of superexchange interactions (in 
10 _3 eV 2 ) in the strong-coupling limit. 



compound 


CF splitting 


TX 


Aai2 


TX 


Aai3 


YTiOs 


-69, -42, 112 


20 





26 


23 


LaTiOs 


-49, 5, 44 


51 


19 


57 


35 



TABLE II: Magnetic interactions (J, in meV) and total ener- 
gies (E, in meV/f.u., measured from the experimentally ob- 
served magnetic state) obtained in Hartree-Fock calculations 
without spin-orbit interaction. 







YTi0 3 






LaTiOs 




phase 


J12 


J13 


E 


J12 


Jis 


E 


F 


2.0 


0.6 





1.2 


0.2 


3.2 


A 


1.6 


-0.2 


0.5 


0.9 


-5.7 


-2.1 


C 


1.4 


0.2 


6.8 


-0.9 


-1.6 


5.9 


G 


1.2 


-1.3 


5.9 


0.1 


-4.6 






robust even in the distorted perovskite compounds. As 
we will see below, it may help to explain the experimen- 
tally observed magnetic ground state in YTO (but not in 
LTO). (iv) The SO interaction at the Ti sites, £~23 meV, 
is also comparable with the CF splitting, and exceeds the 
total energy difference between different magnetic states 
( Table ITT|) . Therefore, it is reasonable to expect essen- 
tially noncollinear magnetic GS with a considerable con- 
tribution of the orbital magnetic moments. 21 

All these trends are clearly seen in HF calculations, in 
which the one-electron TB Hamiltonian for the t<i g bands 
was combined with the on-site Coulomb and exchange 
interactions extracted from the constraint-LDA calcula- 
tions (unless it is specified otherwise)^ and (optionally) 
the SO interaction. The HF potential was treated in the 
rotationally-invariant form™ 

It is true that both in LTO and YTO, the OO is 
strongly constrained by the lattice distortion so that the 
visual change of the OO is not particularly strong in the 
row of FM, A-, C-, and G-type AFM states (Fig. 01. The 
basic question, however, is how this change is reflected 
in the change of other parameters. Our main concern 
is the behavior of inter-atomic magnetic interactions J\y 
Since Jij may depend on the magnetic state (through the 
change of the OO), Eq. (0) is no longer valid. Instead, we 
evaluate Jy separately for different magnetic states using 
the second derivatives of the total energy with respect to 
the angles between spin magnetic moments™ The re- 
sults summarized in Table HT1 clearly show that even tiny 
change of the OO may produce a dramatic change of 
Jij. In addition to the A-type AFM ordering, expected 
from the lattice distortion, the FM state (Ji2>0, Ji3>0) 
can be stabilized by the K&K mechanism both in YTO 
and LTO. 23 Since Ji2>0, the G-type AFM state is un- 
stable. In LTO it can be stabilized only for /7~4.5 eV 
(which leads to Ji2=— 0.3 and Ji3=— 3.4 meV). However, 
this U will also destroy the FM GS in YTO (J 12 =-0.5 
and Ji3=— 0.7 meV). Thus, there is no such parameter 
U which would account for the experimental behavior of 
both YTO and LTO on the level of mean-field HF cal- 
culations. Contrary to the experimental finding^^ the 
magnetic interactions are strongly anisotropic. 

The SO interaction gives rise to a noncollinear mag- 
netic ordering. 21 However, it does not solve the prob- 
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FIG. 3: t2^-electron densities in the ferromagnetic (F), A-, C- 
, and G-type AFM states of LaTi03, without SO interaction. 
Different spin sublattices are shown by different colors. 

lems of the HF description. The magnetic GS real- 
ized in YTO is G a -Ab-F c , which is consistent with 
the neutron-scattering data. 2 Both spin (M5) and or- 
bital (M^) magnetic moments have nonvanishing pro- 
jections onto all three orthorhombic axis a, b, and c, 
which are ordered according with the G-, A-, and F- 



type, respectively. The vectors themselves are given by 
(in fji B , refereed to the site 1): M 5 =(0.05,0. 83,0.34) 
and M L =(-0. 23, -0.33,0.03). The relative weight of 
the F and A components in this structure is very sen- 
sitive to the value of U. The F component will dom- 
inate for smaller U, due to the enlarged K&K contri- 
bution to the 00: e.g. M 5 =(-0.07,-0. 14,0.96) and 
M L =(0.17,0.14 -0.08) for Z7=2.5 eV. The magnetic GS 
obtained in LTO on the level of HF calculations is C a -i^b- 
A c , which has large A component along the c direction: 
M 5 =(-0.13,0.18,0.89) and M L =(-0.14 -0.07 -0.21). 
The G-type AFM structure is totally excluded from C a - 
F\>-A c . Therefore, there is a qualitative inconsistency 
between results of HF calculations and the experimen- 
tal data for LTO. Formally, the problem can be re- 
solved by using larger Z7=4.5 eV, which enforces the 
strong-coupling limit (Table [I]) and leads to the new 
magnetic GS: A a -G h -C c with M 5 =(0.31,0.88,-0.14) and 
M L =(-0. 19, -0.23,0.04). However, the same U would 
lead to the new magnetic GS also in YTO: C a -F\)-A c 
with M 5 =(0.11 -0.19,0.72) and M L =(0.09,0.12,-0.09), 
in disagreement with the experiment. 2 

In summary, the lattice distortion alone does not pro- 
vide a coherent explanation for the unusual magnetic 
properties of YTO and LTO. The complexity of these 
compounds is related with the fact that the CF splitting, 
the SE and SO interaction energies are of the same order 
of magnitude, and should be treated on an equal footing 
beyond the mean-field HF approximation. 
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